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1  echnical  Notes 


Water  Wave  Attenuation  Over  Nearshore 
Underwater  Mudbanks  and  Mud  Berms 

Purpose 

A  shallow-water  wave-mud  interaction  model  is  described  and  used  to 
discuss  design  parameters  for  berms  constructed  from  fine-grained  sedi¬ 
ment.  Model  application  is  presented  to  determine  the  stable  elevation  of 
the  berm  crest  and  the  water  column  height  above  the  crest  in  a  given 
coastal  environment.  The  model  will  predict  wave  attenuation  and  consid¬ 
ers  water  to  be  inviscid  and  mud  to  be  a  highly  viscous  fluid.  It  is  the 
simplest  model  applicable. 
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Background 


Naturally  occurring  underwater  mudbanks  are  known  to  absorb  water 
wave  energy  and  thereby  attenuate  w’aves  that  pass  over  them.  Energy  re¬ 
ductions  on  the  order  of  30  percent  to  as  much  as  90  percent  are  common 
even  in  the  absence  of  any  mea.surable  wave  breaking.  In  recent  years,  en¬ 
gineering  efforts  have  made  use  of  this  property  by  creating  underwater 
mud  berms  to  lessen  wave  impact  in  the  area  Icew’^ard  of  the  berm.  Thus, 
by  appropriately  placing  fine-grained  dredged  material  from  navigation 
channels,  the  disposed  material  can  be  made  to  serve  benefidaliy. 


The  capacity  of  naturally  occurring  bottom  mud  to  absorb  wave  energy 
in  coastal  areas  has  been  known  for  years,  and  in  recent  times  this  prop¬ 
erty  has  been  recognized  as  a  potential  means  for  protGcting  beaches  that 
are  .susceptible  to  erosion  under  w'^ave  attack.  As  an  example  of  the  for- 
nt  mudbanks  occur  on  the  eastern  margin  of  the  Louisiana 
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chenier  plain  where  episodic  waves  are  measurably  damped  as  they  pass 
shoreward  over  the  banks.  The  alongshore  length  of  the  banks  varies 
from  1  to  5  km,  width  from  0.5  to  3  km,  and  thickness  from  0.2  to  1.5  m, 
in  w'ater  depths  from  1  to  2  m  above  the  mud  surface.  The  median  parti¬ 
cle  size  is  3  to  5  pm  (Wells  1983,  Wells  and  Kemp  1986). 


As  an  example  of  a  nonsacrificial,  engineered  berm  off  Dauphin  Island 
in  Alabama,  an  underwater  mound  composed  of  sediment  dredged  from 
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the  Mobile  Ship  Channel  within  Mobile  Bay  was  deposited  by  the  U.S. 
Army  Engineer  District,  Mobile,  in  1988.  The  design  dimensions  of  the 
-berm  were  2,750  m  alongshore,  300  m  wide  at  the  toe,  and  6  m  high  with 
6.5  m  mean  depth  of  water  above  the  berm  crest  (McLellan,  Pope,  and 
Burke  1990).  The  placed  material  is  highly  graded,  with  a  median  size  of 
about  15  |im. 


The  degree  to  which  such  natural  or  manmade  coastal  features  attenu¬ 
ate  incoming  wave  energy  depends  on  their  diiuensions,  the  composition 
of  the  sediment,  and  the  incoming  wave  characteristics.  Thus,  a  predic¬ 
tive  method  for  berm  design  must  include  the  essential  physical  modes  of 
interaction  between  the  flow  and  the  dissipative  bottom  material.  The 
two  most  important  berm  design  parameter  are  tlic  crest  elevation  and 
the  depth  of  water  above  the  cre.st,  since  these  measurably  influence  the 
decree  of  wave  attenuation,  and  hence,  the  wave  impact  in  sheltered  areas. 


Research  performed  under  the  Dredging  Research  Program  (DRP)  Tech¬ 
nical  Area  1,  Work  Unit  "Cohesive  Sediment  Processes,"  has  included  labo¬ 
ratory’’  tests  on  wave  damping  and  fluidization  of  mud,  and  the  develop¬ 
ment  of  two  wave-mud  interaction  models.  The  first  and  simplest  model 
assumes  the  mud  to  be  a  viscous  fluid  and  the  shallow- water  waves  to  be 
linear,  while  the  second  model  considers  the  mud  to  be  viscoelastic  and 
the  waves  to  be  finite  amplitude  (nonlinear).  The  simple  model  is  used  in 
this  technical  note  as  a  guide  for  berm  design  and  to  demonstrate  the  gen¬ 
eral  approach. 


Additional  Information 

This  technical  note  is  based  on  work  by  Drs.  Ashish  Mehta  and  Feng 
Jiang  at  the  University  of  Florida  under  contract  to  the  DRl^.  For  addi¬ 
tional  information,  contact  Mr.  Allen  Teeter,  (601)  634-2820,  or  the  DRP 
manager,  Mr.  E.  Clark  McNair,  (601 )  634-2070. 


General  Physical  Setting  and  Processes 


The  physical  setting  relative  to  the  mudbank  or  the  underwater  berm  is 
conceptually  shown  in  Figure  1,  Seawater  {density  p,  dynamic  viscosity 
}i)  surrounds  the  mud  berm,  winch  is  sn  genera!  a  viscoelastic  material 
(density  p/^j)  possessing  both  an  elastic  component  characterized  by  the 
shear  modulus  G  (typically  on  the  order  of  10^  to  10“^  Pa)  and  viscosity 
(typically  on  the  order  of  lO’  to  Pa  •  sec).  (A  notation  list  is  provided 
as  Table  1.)  Most  muds  are  pseudoplastics  in  terms  of  the  relation  be¬ 
tween  the  shear  stress  and  the  rate  of  strain,  so  that  thev  exhibit  crcco 

A, 

even  at  I'ery  low  applied  stre.s.ses.  Dredging  Research  Technical  Notes  DRP-2- 
04  (Teeter  1992)  describes  the  general  aspects  of  channel-bottom  mud  \ds- 
cous  characteristics. 
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Figure  1.  Schematic  sketch  showing  wave  propagation  over  an  underwater  mud  berm 


When  subjected  to  wave  action,  bottom  muds  respond  by  oscillating  at 
the  forcing  wave  frequency.  As  a  result  of  high  viscosity,  the  oscillations 
attenuate  much  more  rapidly  with  depth  within  mud  than  in  the  water 
column  above.  A  high  rate  of  energy  dissipation  within  the  mud  causes 
the  surface  wave  height  to  decrease  rapidly  with  onshore  distance.  Thus, 
given  the  surface  wave  amplitude  at  .r  =  0  at  the  seaward  edge  of  the 
berm  crest,  the  amplitude  at  any  distance  x,  that  is,  exp  (-kjx),  de¬ 

pends  on  the  wave  attenuation  coefficient  k;.  If  the  bottom  were  rigid,  kj 
would  be  on  the  order  of  10'^  whereas  for  mud,  values  as  high  as  lO"^ 
to  10’^  nT^  are  common  Oiang  and  Mehta  1992).  The  effect  is  much  re¬ 
duced  wave  activity  in  the  area  leeward  of  a  mud  berm  compared  with 
areas  where  the  bottom  Is  composed  of  hard  material.  On  the  other  hand, 
the  wave  celerity  is  higher  where  mud  is  the  substrate  because  the  avail¬ 
able  depth  for  wave  propagation  is  effectively  greater  than  the  water 
depth,  h.  Consequently,  the  depth  h^,  shown  in  Figure  1,  requires 
specification. 

Mud  oscillation  primarily  occurs  as  a  result  of  w'ave-induced  pressure 
work  within  the  body  of  the  material,  while  the  effect  of  shear  stress  is 
more  important  at  the  mud  surface  where  it  can  cause  particulate  resus¬ 
pension.  Thus,  under  continued  wave  action,  the  long-term  equilibrium 
water  depth  above  the  crest  is  that  depth  at  which  the  wave-induced  hy¬ 
draulic  shear  stress  amplitude  is  equal  to  the  erosion  shear  strength  T.; 
of  mud. 

Tn  many  naturally  OL'Curring  environments,  mud  oscillations  under  typi¬ 
cal  fair-weather  wave  conditions  occur  without  much  sediment  resuspen¬ 
sion  and  associated  turbidity  (Jiang  and  Mehta  1992).  Thus,  Tj,,,,  =  can 
be  used  as  a  de.sign  criterion  for  the  berm  crest  elevation  and  the  berm 
can  fulfill  its  role  as  a  wave  attenuator  without  generating  excessive  turbid- 
-  ity  and  dispersion  from  the  site  through  transport  of  eroded  sediment. 

The  berm  slope  from  crest  to  toe  is  determined  by  a  characteristic  yield 
strength,  which  for  design  purposes  has  been  approximated  by  the  upper- 
Bingham  yield  strength  of  the  pseudoplastic  material  (Migniot  1968). 
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Table  1.  Notation  Definitions 


Notation 


Definition 


rfaco  wave  amplitude  at  the  offshore  berm  toe  (ir  =  0 


Surface  wave  amplitude  at  some  inshore  distance  x,  m 

. 

Acceleration  of  gravity,  m/ sec 


Shear  modulus,  I’a 


Water  denth,  m 


Mud  depth  or  tiuckness,  m 


Wave  boundary  layer  thickness  in  mud,  m 
Wave  number,  m'^ 

Wave  attenuation  coefficient,  m'* 

Starting  time,  day 


:  Equilibrium  time,  days 
j  Wave-induced  water  velocity,  m/sec 
Wave-induced  Stokes'  drift  in  mud,  m/scc 
Wave-induced  mud  velocity,  m/sec 
Distance  inshore  of  offshore  berm  toe,  m 


Friction  coefficient.  Pa  ■  sec^/ 

Slope  angle,  deg 

Empirical  slope  constant,  deg/Pa 
Normalized  density  jump  (pM  -  p)/ pM 


:  Nondimensional  wave  number 
Water  surface  deviation,  m 
Mud  surface  deviation,  m 
j  Dynamic  viscosity  of  water.  Pa  •  sec 


Dynamic  viscosity  of  mud.  Pa  •  sec 


Kinematic  viscosity,  m^/ sec 
Nondimensional  wave  attenuation  coefficient 


I  Water  density,  kg/mi 
Mud  density,  kg/m'’ 

Wave  angular  frequency,  rad /sec 


T(,,„  Hydraulic  shear  stress  at  mud  surface,  l^a 

Tc  Erosion  threshold  for  mud.  Pa 

Ty  Upper-Bingham  yield  stress,  Pa 

X  Nondimensional  mud  depth 
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However,  the  stability  of  the  berm  crest  is  not  ensured  solely  through  this 
criterion  for  erosion  since,  due  to  ooen  oarticle  orbits  (Fieure  2)  arisine 
from  nonlinear  wave  effects,  the  residual  velocity  {Stokes'  drift)  can 
cause  the  mud  mass  to  be  transported  landward.  The  impetus  for  tltis  mo¬ 
tion  is  the  net  wave-induced  thrust  that  occurs  in  the  mud  due  to  the 


rapid  wave  attenuation  with  onshore  distance.  Thus,  hypothetically  start¬ 
ing  at  time  f  with  a  mud-water  interfacial  profile  that  can  be  shown  to  be 
near-exponential  in  form  (Jiang  1993),  the  depth-averaged  value  of  Ui  will 
become  nil  everywhere  at  some  instant  tg,  when  an  equilibrium  interfacial 
profile  is  established.  Such  a  condition  w’ill  result  from  a  balance  between 
the  w'ave-induccd  thrust  and  the  adverse  hydrostatic  gradient  in  the  pres¬ 
ence  of  a  sloping  bottom,  as  shown  in  Figure  2.  Profiles  such  as  these  are 


seasonally  established,  for  instance,  along 


the  southwest  coast  of  India  off 


rala,  where  the  southwesterly  monsoon  wav'es  push  the  mudbaTiks 


from  about  20-m  water  depth  to  a  position  close  to 


the  shoreline. 
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where  wave  spectra  were  obtained  at  tw'o  sites  (Mathew  1992),  The  mean 
depth  of  water  at  the  offehore  site  w'as  10  m,  and  at  the  inshore  site  a  lit¬ 
tle  over  5  m.  In  fair  weather,  the  mudbank  off  Kerala  stays  w'^ell  seaward 
of  the  offshore  site. 


Figure  4  is  a  schematic  profile  of  the  Kerala  mudbank  at  its  nearshore  lo¬ 
cation  (Nair  1988).  Its  position  has  been  stabilized  against  the  com¬ 
paratively  steep  nearshore  bottom.  During  two  particular  monsoon  peri¬ 
ods,  the  shoreward  edge  of  the  mudbank  w'as  landward  of  the  inshore 
measurement  site,  while  the  seaw'ard  edge  of  the  mudbank  was  landward 
of  the  offshore  iheasurement  site.  Thus,  the  bottom  at  the  offshore  site 


was  consistently  devoid  of  mud.  Figure  5  shows  examples  of  wave  spec¬ 
tra,  with  and  without  the  mudbank  present,  at  the  inshore  site.  In  fair 
weather,  the  energy  reduction  was  negligible.  (Note  the  slight  phase  shift 
of  the  inshore  spectrum  relative  to  offshore,  which  may  be  an  artifact  of 
the  data  analysis.  If  a  physical  basis  was  the  cause,  it  remains 
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Figure  3.  Coasts!  site  at  AHeppey  in  Kerala,  India,  where  jnor.soonal,  shore-parallel 
mudbanks  occur.  The  pier  is  300  m  long  (after  Mathew  1992) 


Peripheral  Wave 
Damping  Zone 


High  Water  Une 


Figure  4.  Schematic  profile  of  mudbank  region  off  the  coast  of  Kerala,  India 
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Figure  5.  Offshore  and  inshore  wave  spectra  off  Alleppey;  (a)  without  mudbank,  and 
(b)  with  mudbank  present  (after  Mathew  1992) 


in  the  monsoon  case,  the  wave  energy  reduction  v\'as  about  85  percent 
at  the  inshore  site.  In  fact,  as  a  result  of  the  mudbank,  the  wave  energy 
at  the  inshore  site  was  lower  in  monsoon  than  in  fair  weather,  despite  the 
occurrence  of  considerably  more  inclement  offshore  wave  activity  during 
the  monsoon.  Not  surprisingly,  during  monsoon  the  mudbanks  in  Kerala 
have  served  as  safe,  open-coast  havens  for  small  fishing  boats  that  would 
otherwise  require  sheltered  harbors  (Nair  1988). 


The  site  of  the  Mobile  berm, 
in  a  somewhat  idealized  con¬ 
figuration,  is  shown  in  Fig¬ 
ure  6.  Examples  of  offshore/ 
inshore  wave  spectra  are 
given  in  Figure  7  under  two 
offshore  wave  conditions. 
Wave  energy  reductions  were 
significant,  29  and  46  percent, 
respectively,  although  not 
nearly  as  dramatic  as  at 
Kerala  because  the  finer  and 
more  clayey  Kerala  mud  is 
more  dissipative.  Also,  the 
Mobile  berm  is  quite  narrow 
compared  with  the  mudbank. 
Using  a  hydrodynamic  wave 
model,  McLellan,  Pope,  and 
Burke  (1990)  showed  that  an 
assumed  rigid  berm  crest 


rigUfC  6.  CuilStfUCtion  hiie  (cofnciof  for 
dredged  material  placenjent)  for  tbe  Mobile 
berm  (after  McLellan^  Pope^  and  Burke  1990) 
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Figure  7.  Offshore  and  inshore  wave  spectra  at  the  Mobile  berm  site  for  two  different 
wave  conditions  at  the  offshore  site  characterized  by  the  maximum  wave  height, 

(a)  Hniax  =  0-9  and  (b)  H^ax  =  "i-S  m  (after  McLellan,  Pope,  and  Burke  1990) 


produced  negligible  wave  energy  dissipation.  It  is  believed  that  as  a  re¬ 
sult  of  the  depth  of  water,  on  the  order  of  6.5  m  over  the  berm,  the  main 
cause  of  the  observed  damping  is  energy  absorption  by  the  deposited 
mud  at  the  berm  site. 


8 


Technical  Note  URr-1-14  iApJ-i!  1994) 


Interestingly,  diver  observations  at  the  Mobile  berm  site  have  suggested 
the  occurrence  of  surface  wave-forced  interfacial  mud  waves  propagating 
along  the  compliant  crest.  Such  bed  movement  confirms  the  participation 
of  the  bottom  material  in  the  energy  dissipation  process.  Furthermore,  the 
observed  stability  of  the  berm  points  to  the  fact  that  wave  action  becomes 
sufficiently  weak  over  the  berm  to  prevent  significant  wave-induced  scour 
and  degradation  of  the  berm. 

Mass  transport  due  to  Stokes'  drift  is  important  in  situations  such  as 
those  described  earlier  for  Kerala  during  the  Indian  monsoon  when  sus¬ 
tained  and  significant  wave  action  may  occur  for  periods  of  weeks  or 
months.  Furthermore,  the  occurrence  of  a  gravity  slide  into  deeper  waters 
at  the  cessation  of  high  wave  action  toward  the  end  of  the  monsoon 
causes  the  situation  to  be  reinitialized  with  respect  to  mass  transport  at 
onset  of  the  following  monsoon.  Field  evidence  elsewhere,  for  example 
off  Surinam  (Wells  1983),  indicates  that  unlike  the  Kerala  mudbanks  that 
are  clearly  quite  transitory,  mudbanks  typically  tend  to  remain  practically 
stationary  because  of  small  bottom  slopes,  mild  and/or  constant  wave  ac¬ 
tion,  or  other  conditions.  Consequently,  their  net  motion  need  not  be  con¬ 
sidered  in  simple  calculations  for  stable  berm  design. 


A  Simplified  Approach  to  Berm  Design 

Consideration  is  given  to  the  berm  elevation  and  water  depth.  It  is  as¬ 
sumed  that  the  berm  is  endless  in  the  longshore  direction.  This  assump¬ 
tion  requires  the  berm  to  be  much  longer  than  the  length  of  the  wave. 

End  effects  are  therefore  ignored  alongshore,  but  seaward  edge  effects  can 
be  prominent  in  the  cross-shore  direction  if  the  side  slope  is  steep  and  the 
crest  elevation  high  enough  to  cause  an  abrupt  depth  change  at  the  sea¬ 
ward  edge.  Furthermore,  a  change  in  the  sediment  composition  between 
the  ambient  bottom  and  the  berm  may  also  measurably  influence  the 
wave  field.  The  mechanics  of  propagation  of  the  wave  over  the  berm 
would  thus  become  an  initial  (spatial)  value  problem,  the  initial  wave 
amplitude  condition  being  defined  at  the  toe  of  the  berm.  For  typical 
natural  mudbanks,  the  crest  elevations  are  not  high  and  side  slopes  are 
mild.  The  only  real  change  encountered  sometimes  is  the  bottom  composi¬ 
tion,  Even  at  berms  constructed  by  dredged  material  placement,  side 
slopes  are  not  always  very  steep;  at  the  Mobile  site,  the  slopes  range  from 
1:24  to  1:130.  Thus,  for  simplicity,  the  effects  of  both  changing  depth  and 
bottom  composition  on  the  dissipation  process  over  the  berm  can  be  as¬ 
sumed  to  be  negligible. 

Harmonic  solutions  involving  wave  propagation  in  a  single  x  direction 
-  in  water  of  depth  h,  over  mud  of  thickness  hy,  (=  hy),  as  shown  in  Fig¬ 
ure  1,  were  found  for  the  wave  speed  and  other  wave-induced  quantities. 
The  rheological  constitutive  properties  of  mud  are  important,  and  various 
viscoelastic  models  have  been  considered  (such  as  Chou  1989  and  Jiang 
1993).  Recent  studies  on  wave-mud  interaction  by  Feng  et  al.  (1992)  sug¬ 
gest  that  wave  pressure  work  can  cause  the  otherwise  undisturbed  mud 
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to  weaken  through  thixotropic  yield,  leading  to  plastic  behavior  and  loss 
of  effective  stress.  This  process  of  fluidization  is  reflected  in  a  drastic  re- 
■  duction  in  the  shear  modulus  G,  with  the  result  that  the  weakened  mud 
can  be  approximated  as  a  highly  viscous  fluid.  In  this  technical  note,  for 
simplicity,  mud  is  assumed  to  be  a  viscous  fluid  and  water  is  inviscid, 


The  shallow-water  assumption  is  invoked,  and  the  wave  propagation  prob¬ 
lem  »  considered  to  be  linear.  These  two  assumptions  have  been  relaxed 
in  Jiang  (1993)  to  account  for  the  effects  of  finite  amplitude  (nonlinear) 
waves  in  intermediate  water  depths. 


As  a  result  of  the  shallow-water  assumption,  only  horizontal  motion  is 
considered  in  the  governing  equations  of  motion  and  continuity  for  the 
water  layer  (characterized  here  by  no  subscript)  and  for  the  mud  layer 
(subscript  M).  The  horizontal  x  and  the  vertical  y  axes  are  located  at  the 
rigid  bottom;  surface  deviations  q(x,0  and  r\i,Ax,t)  are  functions  of  x  and 
time  t;  u(x,t)  and  uj^(x,y,f)  are  the  wave-induced  velocities;  y  is  the  normal¬ 
ized  density  jump;  hj  is  the  elevation  of  the  water-mud  interface;  and  v  is 
the  kinematic  viscosity  of  mud.  Thus,  y  =  (py^  -  p)/pM/  hj  =  +  t\m^  and 

''m  PM'  where  is  the  dynamic  viscosity. 


Details  of  solutions  for  the  shallow-water  wave  equations  of  motion  and 
continuity  for  the  two-layer  system  including  q,  u,  are  pre¬ 

sented  by  Jiang  (1993).  The  validity  of  the  results  has  been  demonstrated 
by  comparison  to  previous  laboratory  experiments  by  Cade  (1958).  The  re- 
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o  been  applied  to  measurements  obtained  at  a  site  in  the  lit¬ 
toral  margin  of  Lake  Okeechobee,  Florida,  Vv’^here  about  one  third  the  bot¬ 
tom  consists  of  a  mud  layer  of  peaty  origin  (maximum  thickness  0.8  m) 
(Jiang  and  Mehta  1992).  Wave  damping  is  often  significant.  During  a 
field  deployment  it  was  estimated  that  the  wave  heights  were  about  one 
fourth  the  value  that  would  have  occurred  if  the  bottom  had  been  rieid. 


Data  at  the  shallow-water  site  consisted  of  surface  wave,  water  velocity, 
and  mud  acceleration  spectra.  The  water  depth  was  1.43  m,  and  mud 
thickness  0.55  m.  Using  the  wave  spectrum  as  input,  the  velocity  and  ac¬ 
celeration  spectra  were  predicted  and  compared  with  field  data.  The 
model  was  able  to  reproduce  important  features  of  the  velocity  and 
acceleration  spectra.  However,  discrepancies  between  measurements  and 
simulations  occurred  because  of  limitations  in  the  model,  especially  the  as¬ 
sumption  of  a  viscid,  nonelastic  mud  (Jiang  1993),  and  possibly  because  of 
difficulties  in  data  collection  (Jiang  and  Mehta  1 992). 


The  analytic  solutions  for  the  wave  attenuation  coefficient  and  the  wave 
riumber  are  pre.sented  graphically.  In  Figure  8,  the  diniGnsionless  wave  at¬ 
tenuation  coefficient,  ^  where  g  is  gravitational  accelera¬ 

tion  and  a  is  the  angular  frequency  of  the  wave,  is  plotted  as  a  fimction 
of  the  dimen.sion!ess  mud  depth,  x  =  for  ranging 

from  0.1  to  1,  This  solution  is  valid  only  for  y  =  -  p)./p,\^  =  0.15, 

which  is,  however,  a  typical  value.  It  can  be  shown  that  (2py^/p^»^.{CrW"  is 
representative  of  the  thickness  of  the  wave  boundary  layer  in  mud.  For  a 
given  water  depth  and  angular  wave  frequency  o,  attenuation  is  observed 
to  be  maximum  for  all  hj^/h,  when  x  is  equal  or  close  to  unity. 
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Figure  8.  Relationship  between  \  and  %  representing  the  effects  of  wave  energy 
dissipation  in  mud  on  wave  attenuation  in  shallow'  wafer.  (Corresponding  wave 

number  presented  in  Figure  9> 


Therefore,  this  condition  of  maxiinum  kj  means  that  the  rate  of  wave  en¬ 
ergy  dissipation  in  the  mud  is  maximum  when  the  wave  boundary  layer 
attains  about  the  same  thickness  as  the  mud. 

Two  solution  limits  must  be  mentioned.  One  is  x  =  0^  which  corre¬ 
sponds  to  a  mud  of  infinitely  high  viscosity  corresponding  to  a  rigid  bot¬ 
tom.  In  this  case,  there  is  no  participation  by  mud  in  the  dissipation  pro¬ 
cess  and,  since  water  is  assumed  to  be  inviscid,  tliere  is  no  damping.  The 
other  limit  is  approached  with  increasing  x  corresponding  to  decreasing 
mud  viscosity.  In  this  case,  damping  also  approaches  zero  as  suggested 
by  the  display. 

In  Figure  9,  the  companion  plot  for  determining  the  wave  number  k  is 
given.  The  dimensionless  wave  number  ^  =  k/[a/{gh)^^^\  is  shown  as  a 
function  of  %.  At  x  =  0,  ^  =  1  corresponds  to  shallow-water  wave  motion 
in  depth  ft,  whereas  with  increasing  x  mud  depth  increases,  and  in  the 
limit  the  problem  becomes  one  of  inviscid,  shallow-water  w'ave  propaga¬ 
tion  in  depth  ft  +  ft/^.  Note  that  the  solution  for  the  maximum  crest  water 
velocity  amplitude  at  the  .seaward  edge  of  the  berm  crest  is  =  (i{ygk/<y. 
Thus,  the  shear  stress  amplitude  at  the  mud  surface  is  =  a  Um  where 
(x  incorporates  water  density  and  a  drag  term  that  depends  on  the  rough¬ 
ness  of  the  mud  surface  as  well  as  the  amplitude  of  the  water  particle  just 
above  the  mud  surface  (Dyer  1986). 
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Figure  9.  Relationship  between  t  and  x,  for  determining  the  wave  number 

corresponding  to  Figure  8 


Example  Application 


These  solutions  can  be  used  to  address  a  simple,  hypothetical  berm  de¬ 
sign  problem.  Assume  a  3-m-high  ihn/  =  h^)  fine-grained,  dredged  mate¬ 
rial  (density  1,180  kg/m^)  mound  structure  in  water  such  that  it  remains 
stable.  At  the  same  time,  the  economics  of  dredged  material  transport  dic¬ 
tates  that  the  structure  be  placed  in  the  shallowest  possible  depth  to  mini¬ 
mize  the  distance  between  the  shoreline  and  the  structure.  Other  relevant 


parameters  are  wave  frequency  o  =  0.2k  rad/sec,  forcing  wave  amplitude 
(Iq  =  0.75  m,  erosion  shear  stress  =  1  Fa,  and  a  =  1 .03  Pa  •  sec^/m^  for  a 
hydraulically  smooth  turbulent  flow.  First  examine  the  effect  of  mud  vis¬ 
cosity  on  the  problem  by  selecting  two  values  of  the  mud  kinematic  viscos- 
-  Hm/Pm  =  ^  ^0'*  m^/sec.  Then,  find  the  water  depth  h 

and  calculate  wave  attenuation,  assuming  the  berm  to  be  300  m  wide  (lOp 
in  Figure  1). 


For  this  case  the  dimensionless  mud  depth  x  =  53  and  5.3  for  the  low- 
and  high- viscosity  values.  Figures  8  and  9  can  then  be  used  to  find  the 
dependence  of  the  wave  number  k,  the  wave  attenuation  coefficient  k^,  and 
the  bed  shear  stress  amplitude  at  the  berm  crest  X;,,,,  on  the  water  depth  h. 
In  Figure  10,  these  relationships  are  shown  graphically  for  depth  /; 
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Figure  10*  Effect  of  mud  viscosity  on  relationships  relevant  to  berm  design:  (a)  wave 
number  versus  water  depth,  (b)  wave  attenuation  coefficient  versus  water  depth,  and 
(c)  bed  shear  stress  amplitude  versus  water  depth 
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ranging  from  1  to  8  m.  Note  that  for  the  chosen  wave  frequency,  the 
model  cannot  be  used  for  greater  depths  because  of  the  shaiiow-water  as¬ 
sumption.  The  wave  number  is  seen  to  be  comparatively  insensitive  to 
viscosity,  since  4t  depends  primarily  on  the  inertia  and  pressure  forces  in 
the  water  colunin.  The  same  is  observed  in  the  case  of  the  bed  shear 
stress. 


Note  also  that  the  magnitude  of  the  bed  shear  stress  is  an  artifact  of  its 
chosen  dependence  on  the  water  velocity  in  the  inviscid  water  column 
under  the  given  assumptions.  This  dependence  is  characterized  by  the  co¬ 
efficient  cx.  On  the  other  hand,  wave  damping  is  seen  to  be  highly  sensi¬ 
tive  to  the  viscosity,  which  emphasizes  the  need  to  appropriately  consider 
the  composition  and  the  constitutive  behavior  of  the  particular  mud  under 
investigation. 


Given  the  erosion  shear  stress  tj  =  1  Pa,  the  minimum  water  depth  re¬ 
quired  for  a  stable  berm  would  be  3.2  m  (Figure  10c),  using  the  higher  of 
the  two  viscosities  as  most  realistic.  Tlie  site  where  this  depth  occurs  will 
be  the  minimum  distance  from  the  shoreline  where  the  berm  must  be 
placed  for  it  to  be  nonsacrificial.  Since  for  the  depth  of  3.2  m  the  wave 
attenuation  coefficient  ki  -  0,002  1/m,  tiie  1.5-m-high  wave  will  be  re¬ 
duced  to  0.82  m,  which  corresponds  to  a  70-percent  reduction  in  the  wave 
energy. 


Information  on  the  borm  crest  elevation  and  the 
for  a  stable  berm  must  be  supplemented  by  other  design  parameters  foi 
engineering  implementation  of  berm  design.  Among  these,  berm  slope  is 
perb.aps  most  important.  In  that  context  the  relationship  of  Migniot 
(1968),  tan  cc'  =  where  of  is  the  slope  (angle),  is  worth  a  note.  From 
laboratorv  experiments  the  empirical  constant  B  was  found  to  vary  from 

X  X  •  ✓ 

0.007  to  0.025  when  the  yield  stress  Zy  was  measured  in  pascals.  If  for  the 
sake  of  argument  it  is  assumed  that  Zy  =  t,.  =  1  Pa,  the  corresponding 
range  of  slope  will  be  1:143  to  1:40,  which  is  realistic  when  compared  with 
1:130  to  1:24  at  the  Mobile  berm. 


Conclusions 


The  Mobile  berm  appears  to  be  effective  in  damping  approaching  water 
waves  because  of  the  highly  dissipative  nature  of  tlie  mud  bottom.  In 
that  context,  the  model  presented  here  is  a  useful  tool  for  examining  the 
essential  physics  underlying  the  dissipative  mechanism,  which  is  almost 
entirely  embodied  in  the  rheology  of  the  mud.  An  application  of  the 
shallow-water  model  to  calculate  design  berm  crest  elevation  and  water 
depth  w’^as  presented.  The  model  relies  on  a  minimum  of  input  informa¬ 
tion  and  is  amenable  to  graphical  presentation  of  berm  design  calcula¬ 
tions.  In  some  ca.ses,  however,  the  more  complex  model  presented  by 
Jimig  (1993)  may  be  more  appropriate.  It  better  incorporates  physics  of 
the  wave-mud  interaction  problem,  and  can  be  potentially  used  for  berm 
design  calculations  under  conditions  for  which  the  assumptions 


Techtncal  Nofr  DRP-l-14  (April  1994) 


underlying  the  shallow-water  model  are  not  valid.  That  model  can  also 
be  used  to  assess  potential  wave  impacts  where  a  preexisting  berm  or  a 
bank  is  removed  naturally  or  otherwise,  thus  exposing  a  hard  bottom  over 
which  wave  dissipation  would  be  reduced. 
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